Long-range quantum gates using dipolar crystals 
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We propose the use of dipolar spin chains to enable long-range quantum logic between distant 
qubits. In our approach, an effective interaction between remote qubits is achieved by adiabatically 
following the ground state of the dipolar chain across the paramagnet to crystal phase transition. 
We demonstrate that the proposed quantum gate is particularly robust against disorder and derive 
scaling relations, showing that high-fidelity qubit coupling is possible in the presence of realistic 
imperfections. Possible experimental implementations in systems ranging from ultracold Rydberg 
atoms to arrays of Nitrogen- Vacancy defect centers in diamond are discussed. 

PACS numbers: 03.67.-a, 05.30.Rt, 76.30.Mi, 32.80.Ee 



The ability to carry out quantum gates between spa- 
tially remote qubits forms a crucial component of quan- 
tum information processing [l|. Theoretical and exper- 
imental work addressing this challenge has largely been 
focused upon using photons 043], spin chains [B-lal an d 
other hybrid systems 0- 11 1 as quantum buses, which me- 
diate long-range quantum information transfer. In these 
approaches, this transfer is achieved by either encoding 
the information in a traveling wavepacket 043 > or by 
coupling the remote qubits to a shared spatially extended 
mode [8l4llj. In this Letter, we propose a novel approach 
to this outstanding problem and demonstrate that adi- 
abatic driving of a dipolar spin system across a quan- 
tum phase transition can be used to implement a robust 
controlled-phase (CP) gate. 

Our approach is applicable to dipolar spin systems 
1 2l-ll 7T| . composed, for example, of ultracold atoms and 
molecules, or solid-state spin ensembles, where natural 
imperfections invariably lead to disorder. E.g., for spin 
qubits associated with Nitrogen- Vacancy (NV) centers 
in diamond [3 [l9|], the need for both long-range and 
disorder-robust quantum gates is especially evident. De- 
spite room temperature coherence times of ~ 10 ms, the 
weakness of magnetic dipolc-dipolc interactions limits 
NV spacing to ~ 10 nm for effective two-qubit gates. 
Even recently demonstrated sub-wavelength techniques 
[iol ] cannot address individual NV qubits at such small 
separations. Moreover, any solid-state quantum bus de- 
signed to mediate longer ranged interactions will suffer 
from positional disorder due to the difficulty of precise 
nanoscalc implantation. 

Here, we explore a possible solution to creating a quan- 
tum bus within a disordered system. The key element un- 
derlying our proposal is a phenomenon discussed in the 
context of Rydberg atoms as the blockade effect 2lM24 j . 
The simultaneous driving of two excitations within length 
scales shorter than the blockade radius is forbidden as 
strong interactions shift the doubly-excited state away 
from resonance. Hence, within the blockade radius, the 
underlying spatial distribution of the sites is largely ir- 



relevant and the arising many-body ground state washes 
over the effects of disorder and can lead to the formation 
of crystalline structures [l]| . 



The dynamic crystal formation [13Ml6| underlying our 
protocol is schematically illustrated in Fig. 1. We con- 
sider two qubits A and B, coupled to the ends of a one- 
dimensional (ID) quantum bus of two state atoms with 
one electronic ground state and one Rydberg state. The 
quantum bus is initially prepared in its ground state, 
containing no Rydberg excitations. Then, the bus is adi- 
abatically driven into the crystal regime. The resulting 
many-body state has an energy which depends on the 
boundary conditions set by the state of the qubits. Intu- 
itively, this dependence results from a compression of the 
crystal, and hence, a decrease in the distance between two 
Rydberg excitations, or, when the boundary qubits are 
not excited. Under free evolution, this energy difference 
is translated into a phase difference, which entangles the 
qubits. After reversing the adiabatic step, the quantum 
bus returns to its initial state while the qubits remain 
entangled. 

To be specific, we consider an ensemble of strongly in- 
teracting two-state systems described by the Hamiltonian 
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where A represents the detuning, f2 is the Rabi frequency 
and p — 3 for dipolar interactions or p = 6 for van der 
Waals interactions. The interaction strength is charac- 
terized by the coefficient C p and involve projectors onto 

one of the states, P} = ||Xtl = (1 + of)/ 2 - As w <3 wm 
discuss further below, the same Hamiltonian also applies 
to NV centers under appropriate driving. Here, the elec- 
tronic ground state is a spin triplet; thus, the m s = 
state corresponds to the atomic ground state, while the 
m s = 1 state, which possesses a magnetic dipole moment, 
corresponds to the excited Rydberg state (l8. 19l. 2ol 26|. 



An analogous Hamiltonian governs the interactions be- 
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FIG. 1: Setup for the proposed gate, (a) Depending on the 
state of the qubits A and B, the quantum bus (atoms shown 
as black dots) in the crystalline phase will possess a differ- 
ent ground state, where the distance between two excitations 
an is altered. This intuitively corresponds to a boundary- 
condition-dependent compression of the crystal (a' R < an). 
(b) Ground state and first excited state during the control 
sequence. Initially, the quantum bus is prepared in the para- 
magnet (1). Then, the control parameter A is increased adia- 
batically (2), driving the system across the phase transition. 
Once in the crystal the system evolves freely and picks up a 
phase shift depending on the qubit states (3). After the adia- 
batic process is reversed (4), the quantum bus is disentangled 
from the qubits (5). The inset shows the phase diagram of 
the system, with the dashed line depicting the control profile 
dependent on the Rabi frequency fl(t) and the detuning A(t). 



tween the boundary qubits (A,B) and the quantum bus, 
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Since the interaction conserves a A and a%, any entan- 
gling operation between the qubits will be in the form of 
a CP gate. We may drive the quantum bus independently 
of the qubits by ensuring that the resonant splitting of 
the qubits differs from that of the mediating bus; possible 
implementations will be discussed later. 

To derive general scaling properties, we now consider 



a one-dimensional system containing N two-state spins 
within length L. The CP gate protocol consists of an 
adiabatic ramp from the classical ground state into the 
crystal regime for time to, hold for phase accumulation 
for time t v and reverse ramp (another to), resulting in 
a total gate time t g = 2to + t n . There are three factors 
which influence the asymptotic scaling of the fidelity with 
system size: (i) the strength of the effective interaction 
E- m t at the hold point, (ii) the minimum energy gap A g 
(across ramp and qubit sectors) protecting the adiabatic 
evolution, and (iii) the strength of external decoherence 
mechanisms. The interaction energy between the qubits, 
which governs the timcscale of entanglement generation, 
is 
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where E a p refers to the energy of the many-body state 
with the qubits in state \ce) A \(3) B (Fig. [1]). Within the 
continuum limit of a classical crystal, E- mt ~ d 2 /L for 
d -C cifj, where d is the distance between the qubits and 
the ends of the quantum bus. Owing to quantum fluc- 
tuations, the classical crystal cannot be the true ground 
state, and the system is rather described in terms of a 
Luttinger liquid 



27 



|28| . However, such corrections are 
important only in the limit of very large system sizes 
(N ~ 10 9 for typical parameters), where they lead to an 
algebraic decay of the correlation functions 27 [ . 

To analyze the effects of a finite gap and decoherence, 
we consider the contributions of each to the overall error 
of the controlled phase gate, assuming they occur inde- 
pendently [Isj]. While, in the thermodynamic limit, the 
gap vanishes at the phase transition, here, we consider 
finite system sizes where there always exists a non-zero 
gap. For gapless phases such as the dipolar crystal how- 
ever, it is important to note that the gap may further de- 
crease upon entering the ordered phase. The qualitative 
effect of such a finite gap is described within a Landau- 
Zener framework. Optimizing the form of the Landau- 
Zener sweep by introducing a nonlinearity results in an 
improved error scaling with, elz = exp ( — c A n tp/h) , 
where c is a mo del- dependent numerical constant 
In the case of the dipolar crystal, Ag ~ 1/L, due to the 
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phononic nature of the excitations [2j 
in an error which scales according to the theoretical op- 
timum given by the Lieb-Robinson bound for the speed 
of information transfer [32^ . 

Next, we consider the effects of decoherence, noting 
that the induced error is a monotonically increasing func- 
tion dependent only on the product of the decoherence 
rate 7 and the total gate time t g . In particular, we as- 



sume, Ed = 1— cxp[— (7*3) ] ~ i'ytg) 6 , where S depends on 
the physical details of the decoherence process [33| . The 
highly entangled nature of the various many-body states 
depicted in Fig. Q] implies that the effective decoherence 
rate must scale with the system size, 7 = loj^i where 70 
is the single particle decoherence rate and the length scale 



3 



Lq. In the dipolar crystal, L$ is approximately given by 
the average distance between two excited spins; this is 
consistent with intuition, as decoherence processes ought 
only be relevant at sites where there exists an actual ex- 
citation. 

To separate off the explicit system size dependence 
within Slz, we define ao = cLAq/H. Combining the two 
error contributions and conservatively plugging t g > to 
into the form for clz, then yields 
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By minimizing the total error, we obtain an optimal gate 
time, t° pt = SL\og[L ao/(L 2 j )]/a , with corresponding 
error, 
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Thus, our protocol exhibits a scaling analogous to a quan- 
tum gate based on a microscopic interaction with energy 
C2/L 2 , which has an error given by e = L 2S / C2Y ■ 
The precise values of ao and Lq can be derived from nu- 
merical studies, as we will show in the following. 

Simulations. — In these numerical simulations, we both 
verify our general scaling arguments and demonstrate the 
ability to achieve superior gate fidelities in comparison 
to bare microscopic interactions. We begin by consid- 
ering a chain of N equidistant or uniformly randomly 
distributed particles with average intcrparticle spacing 
a; the details of the numerical simulation method are de- 
scribed in fl3| . Initially, the qubits are prepared in the 
(It) a + I4->a) ® (lt> fl + ll) s )/2, while the 
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spin chain is fully polarized, \^) sc = njt)i- This initial 
qubit state consitutes a worst-case scenario leading to a 
minimum value of the fidelity for relevant decoherence 
models, including pure dephasing. Using a different ini- 
tial state for the qubits can lead to a significantly higher 
fidelity. To drive the system across its phase transition, 
external control fields arc then varied according to 



Q(t) = fin sin 2 , 

\l + 16t 2 /4 

A(t) = A [l-5exp(-4t/t ) 
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(6) 
(7) 



Here, fin is chosen such that the endpoint of the ramp lies 
just within the crystalline phase, see Fig. [T] While the 
proposed ramp profile features the requisite nonlinear- 
ity, its details have yet to be optimized; therefore, with 
optimal control theory, it may be possible to further en- 
hance the achievable gate fidelities j34|. At t = to, the 
system freely evolves for a time, t n = nh/Ei n t, in or- 
der to allow the effective interaction to generate a phase 
gate between the qubits. Following this period of free 
evolution, the adiabatic ramp is then reversed. In addi- 
tion to naturally following the reversed profile, an alter- 
nate implementation can also be achieved by a complete 
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FIG. 2: Numerical simulation results for the gate fidelity Flz 
due to non-adiabatic transitions depending on the product of 
the gap Ac and duration of the adiabatic step to for equidis- 
tant (N — 34, diamonds) and disordered (N = 15, crosses) 
configurations. The gap has been varied by changing fin while 
all other parameters are held fixed. The solid line is an expo- 
nential fit to the data, (p = 3, C3 = lOOMlna 3 , Ao = 2.3 fio, 
d = 3 a). 



reversal of Hamiltonian dynamics upon switching H to 
— H . The ability to change the sign of the interaction 
depends on the physical implementation; in the case of 
Rydberg atoms, this can be achieved by transferring the 
population from a repulsive state to an attractive state 
or by changing an external electric field 



351. We find 



that both protocols give nearly identical results and fo- 
cus on the latter, as it simplifies the numerical analysis. 
Since our procedure implements a controlled phase gate 
up to local rotations, the fidelity of the proposed gate is 
then given by the disentanglement fidelity between the 
qubits and the chain, Flz = yTr {Pab}> where pab is 
the qubits' reduced density matrix. 

In the absence of decoherence, we expect the er- 
rors to be characterized by a Landau-Zener exponen- 
tial and this is indeed revealed by the simulations, as 
shown in Fig. 2; the fidelity is well-characterized by 
F = 1 — 6cxp(— cAQto/h), where b and c are numeri- 
cal fit parameters. Combining these numerics with the 
additional errors associated with decoherence (ed) yields 
an overall fidelity, 
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Note that for near-perfect gates this expression is equiv- 
alent to Eq. (jl]). As previously discussed, there now ex- 
ists a maximum fidelity, which is achieved by an optimal 
ramp time, to that is a function of only the effective in- 
teraction strength, E int , and the gap, Aq. 

Crucially, Eq. (|SJ) now allows us to investigate the con- 
sequences of a disordered interparticle spacing. For a 
ID dipolar crystal, it is known that the crystal spacing, 
a R = [C(p)(p + l)Cp/A] 1 ^ p , is essentially independent of 
the spacing between individual particles, suggesting that 
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FIG. 3: Dependence of the maximum fidelity of the proposed 
quantum gate on the single particle decoherence rate 70 with 
parameters taken from the numerical simulation for a system 
size of L and a spin-echo suppressed decoherence rate with 
S ~ 3 . The solid red line is the fidelity in the equidistant case 
with TV = 34 particles, while the shaded areas correspond to 
90% confidence intervals for a disordered situation (N = 51). 
The dashed line indicates the fidelity that can be achieved 
using the bare dipolar interaction between the qubits. 



the crystalline phase may be robust against effects of dis- 
order 27). To evalutate the gate fidelity (|8]), we numeri- 



cally determine the gap, Ag, and the interaction energy, 
E- m t, for 100 different uniformly distributed random con- 
figurations for the parameters as in Fig. [5] (b = 0.62, c = 
0.32). By extracting L = = 2.0 (C 3 /^ ) 1/3 

from the average density of excitations and employing 
Eq. l[5|). we then calculate the optimum adiabatic ramp 
time to and hence, the maximum gate fidelity. The re- 
sult is shown in Fig. [3] Significantly, even in the presence 
of disorder, the fidelity of our phase gate is higher than 
that which can be achieved via microscopic dipolar inter- 
actions. 

Experimental Realization. — In the case of Rydberg 
atoms, it is possible to achieve independent addressing of 
the boundary qubits by choosing differing hyperfme lev- 
els for the qubit and quantum bus atoms. Our proposed 
protocol can be realized through either van der Waals 
interactions (jp = 6) between S states or through dipolar 
interactions (p = 3) between states within an electric- 
field induced Stark fan [3r| . Let us focus on the latter 
case and consider a Rydberg state with principle quan- 
tum number n = 43 and decoherence rate 70 = 10 kHz 
(S = 1), neglecting the effects of atomic motion, which 
is well justified at temperatures on the order of 100 nK. 
In order to achieve a gate fidelity, F = 0.9, the requi- 
site laser parameters are given by f^o = 2tt x 8 MHz, 
Ao = 27t x 17 MHz, and the corresponding interaction 
strength is given by C?,a~ 3 — 2tt x 800 MHz. Such an 
interaction strength can be achieved in a Rydberg atom 
cloud with an average intcrparticlc spacing of a ~ 1/xm, 
leading to L w 40 /im; we note that these parameters 
arc compatible with present experimental techniques [37| ■ 



For van der Waals interactions, the enhancement of the 
phase gate fidelity is even more pronounced. 

It is also possible to implement our protocol using a 
solid-state room temperature setup based upon NV de- 
fect centers in diamond 0, [l9| . To realize the Hamil- 
tonian (TTJ) using a quasi- ID chain of NV centers (38l - l40| . 
we work with the m s = 0,1 electronic spin states with 
crystal field splitting « 2.8 GHz. Independent addressing 
of the qubits and quantum bus can be achieved by using 
a different Nitrogen isotope ( 15 N) for the qubits and the 
bus ( 14 N); the hyperfme coupling between the NV elec- 
tronic spin and the Nitrogen nuclear spin is isotope de- 
pendent, with A^ 15 « 3.03 MHz and w -2.14 MHz 
4l| . This difference ensures that the microwave driving 



of the quantum bus is off-resonant with the qubit split- 
ting, allowing for independent initialization and control. 
Alternatively, the boundary qubits may be taken as bare 
Nitrogen impurities manipulated by a nearby NV center. 
The Nitrogen electron spin functions as the boundary 
spin for the protocol and features resonance frequencies 
detuned by GHz, which potentially allows for stronger 
driving during the gate procedure. 

While recent experiments have demonstrated optical 
initialization of the m s = state with approximately 92- 
95% fidelity, it may be possible to further enhance this 
initialization by exploiting the neutral NV° charge state. 
In particular, recent work 42j has shown that red laser 
excitation can transfer nearly 100% of the population to 
one m s sublevel of the spin- 1/2 NV° electronic spin. This 
enables us to effectively polarize the NV nuclear spins of 
the chain via cycles of microwave and red laser driving. 
Finally, after returning the defect to the NV~ charge 
state, a SWAP gate transfers the polarization from the 
nuclear spins back to the spin-1 electronic spins of the 
NV~. Such enhanced initialization may prove benefi- 
cial for other NV-based quantum computing architec- 
tures m. 

Magnetic field fluctuations (e.g., from a nuclear spin 
bath), which give rise to T 2 * dephasing processes can ef- 
fectively be suppressed by stroboscopically switching the 
system between the m s = ±1 states. In the resulting dy- 
namics, the electronic spin is then decoupled from the en- 
vironment and coherence times up to the spin relaxation 
time Xi can be achieved . This procedure also leads to 
the suppression of undesired flip-flop couplings between 
the NV centers. Assuming 70 = 100 Hz [44j, a Rabi fre- 
quency Slo = 27r x 62 kHz, a detuning Ao = 2-7T x 130 kHz, 
and an average NV spacing of a — 2 nm, according to 
the results shown in Fig. |3l we can achieve gate times 
t g £3 500 /is and fidelities of F = 0.98 in the equidistant 
case, and F = 0.93 ± 0.04 for disordered configurations 
over a distance of L = 74 nm. We stress that such qubit 
distances are compatible with the individual optical ad- 
dressing and readout of NV qubits using experimentally 



demonstrated sub- wavelength techniques [20l |45|. 



In summary, we have shown that a robust long-range 



5 



quantum gate can be created using dipolar spin chains. 
We have discussed possible experimental realizations 
with Rydberg atoms or NV centers and emphasize that 
the proposed long-range gate can tolerate disorder. At 
the same time, the proposed gate is not limited to the 
case of dipolar crystals; indeed, one can implement our 
protocol within the transverse Ising model, wherein even 
a nearest-neighbor interaction can be used to create an ef- 
fective 1/L 2 power law interaction. Finally, our proposal 
is also a first step towards studying quantum many-body 
physics with NV centers. 
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